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We propose a guideline for exploring substrates that stabilize the monolayer honeycomb structure 
of silicene and germanene while simultaneously preserve the Dirac states: in addition to have a 
strong binding energy to the monolayer, a suitable substrate should be a large-gap semiconductor 
with a proper workfunction such that the Dirac point lies in the gap and far from the substrate 
states when their bands align. We illustrate our idea by performing first-principles calculations for 
silicene and germanene on the Al-terminated (0001) surface of AI 2 O 3 . The overlaid monolayers 
on Al-terminated Al 203 ( 0001 ) retain the main structural profile of the low-buckled honeycomb 
structure via a binding energy comparable to the one between silicene and Ag(lll). Unfolded band 
structure derived from the /c-projection method reveals that gapped Dirac cone is formed at the 
K point due to the structural distortion and the interaction with the substrate. The gaps of 0.4 
eV and 0.3 eV respectively for the supported silicene and germanene suggest that they may have 
potential applications in nanoelectronics. 

PACS numbers: 71.20.-b,73.20.-r,73.22.Pr 


Silicene, germanene and stanene have attracted con¬ 
siderable attention in recent years since the discovery of 
graphene for they can host a variety of novel physical phe¬ 
nomena, such as quantum spin Hall effecfP ^, qu antum 
anomalous Hall effecfP^^ and valley Hall effectP^ol, which 
are related to the Dirac electrons in these unique two- 
dimensional systems. Moreover, t hey po tentially exhibit 
extraordinary electronic properties^^^^, holding signifi¬ 
cant promise for nanoelectronic devices. 

However, in sharp contrast to graphene which can 
be mechanically exfoliated from graphite silicene, ger¬ 
manene and stanene have to be grown on substrates. Up 
to now, silicene monolayers and multi-layers have been 
obtained by molecular-beam epitaxial (MBE) growth on 
Ag(lll)|2J^^, Ir(lll)p2, and ZrB 2 ( 0001 )P^. These sys¬ 
tems exhibit a variety of structural reconstructions such 
as (4x4)P^, (\/l2 X \/l2), (\/l3 x \/l3) with respect to 
the unit cell of Ag(lll) and (\/3 x \/3)P^ with respect 
to the 1x1 silicene. Most recently, germanene^^HMl 
stanen^^ were also obtained by MBE growth. However, 
the structural reconstruction and the strong hybridiza¬ 
tion between the grown monolayers and the substrates 
may sign ificantly modify the electronic structures of the 
overlayer^^SEoJUl] Although theoretical calculations indi¬ 
cate that the Dirac states in the monolayers can be well 
preserved when t hey in teract with the substrate via vdW- 
type interactions^^E^H^, the weak layer interactions are 
unfavorable for stabilizing the low-buckled honeycomb 
structure. On the other band, it is predicted that the 
monolayer structure can be stabilized on semiconducting 
surfaces such as the ( 111 ) surface of GaP, GaAs, ZnS, 
and ZnS^^. However, electronic bands lose features 
of the Dirac states. Therefore, it is essentially important 
to explore substrates that stabilize the monolayer honey¬ 
comb structure while simultaneously preserve the Dirac 
states in the silicene/germanene monolayer. 

we propose that a suitable substrate for silicene and 


germanene should satisfy the following criteria: (i) it has 
a strong binding energy to the monolayer; (ii) it has a 
large energy gap; (iii) it has a proper workfunction such 
that the Dirac point of the monolayer lies in the gap and 
far from the substrate states when their bands align. We 
illustrate our idea by performing first-principles calcula¬ 
tions for silicene and germanene on the Al-terminated 
( 0001 ) surface of AI 2 O 3 . Despite the binding energy be¬ 
tween the monolayer structure and AI 2 O 3 comparable 
to that for metal-supported systems, unfolded bands re¬ 
veal that the Dirac states in silicene and germanene are 
basically preserved on Al-terminated AI 2 03 ( 0001 ). The 
structural distortion along with the interaction from the 
substrate induce a gap opening at K, which makes them 
promising for potential applications in nanoelectronic de¬ 
vices. Moreover, our results reveal that the interaction 
from the substrate gives rise to minigaps in the band 
structure of the overlaid monolayers, which will lead to 
Eermi velocity renormalization and thus affect the elec¬ 
tronic transistor properties. 

Our concept is illustrated in Eig. When a semicon¬ 
ductor substrate and a monolayer with Dirac states are 
brought together, their bands are aligned to a common 
vacuum level. The interaction of between them is de¬ 
noted by Yint- Similar to the heteronuclear diatomic sys¬ 
tems the resulting perturbations to the electronic bands 
of the isolated systems are determined by both the in¬ 
teraction and the difference (AE) in the energy levels 
between the two constituents. Strong perturbations to 
the Dirac states can be expected when Yint is large and 
AE is small. In particular, the silicene/germanene bands 
are strongly destroyed when it is placed onto metal sub¬ 
strates because of the strong layer interaction and that 
AE is zero. Whereas, the Dirac states can be preserved 
well in vdW heterostructures since Yint is small. How¬ 
ever, to stabilize the monolayer honeycomb structure a 
strong binding between silicene (germanene) and the sub- 
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FIG. 1. Schematic illustration of the interaction between a 
semiconducting substrate and a monolayer with Dirac states. 
Wsnb (Wmono) represents the workfunction of the substrate 
(the monolayer). AE^ denote the energy differences between 
the substrate states and the Dirac point. Yint stands for the 
interaction between the substrate and the monolayer. VBM 
and CBM represent the valence band maximum and conduc¬ 
tion band minimum, respectively. 



FIG. 2. Structure of (\/l3 x \/T3) silicene on (3x3) Al- 
terminated AD 03(0001). (a) Perspective view and (b) top 
view of the relaxed SI structure. Top view of the relaxed 
(c) S2 and (d) S3 structures. Si atoms are represented by 
blue, green (buckling away from the substrate), and yellow 
(buckling towards) balls. 


strate is required, that is, Yint is large. Therefore, to 
avoid a strong disturbance to the Dirac states AE should 
be large, that is, the Dirac states should be far from the 
substrate states. For such a purpose, large-gap semicon¬ 
ductors with a proper workfunction are desirable such 
that the Dirac states lie in the middle of the gap when 
their bands align. Of course, the interaction with the 
substrate may affect the structure and correspondingly 
modifies the electronic bands of the overlayer. 

As a large-gap semiconductor, AI 2 O 3 has been widely 
used as a substrate in silicon-based semiconducting de¬ 
vices. A simple density-functional theory (DFT) calcula¬ 
tion indicates that Al-terminated AD 03 ( 0001 ) is a semi¬ 
conductor with a gap about 4.84 eV and a work function 
of about 6.60 eV. While the ideal low-buckled silicene 
has a workfunction of about 4.76 eV. By the band align¬ 
ment, one can find that the Dirac cone of the silicene 
lies in the gap of Al-terminated AD 03 ( 0001 ), and is 1.94 
ev and 2.90 eV from the VBM and CBM, respectively. 
Such energy separations imply that the Dirac states in 
silicene may experience a relatively weak perturbation 
(compared to silicene/Ag(lll)) even though the mono- 
layer has a strong binding with the substrate as well. 

To validate our ideal, we performed DFT calculations 
for silicene on Ap 03 ( 00 01 ) using the Vienna Ab initio 
Simulation Packagd^®^. The pseudopotentials wer e con- 
structed by the projector augmented wave methocP^EH. 
Van der Waals dispersion forces between the adsor¬ 
bate and the substrate were accounted for through the 
optPBE-vdW functional by using the v dW-D F method 
developed by Klimes and Michaelide d^^ l ^^ i A 4x4 
Monkhorst-Pack k-mesh was used to sample the surface 
BZ and a plane-wave energy cut off of 400 eV was used 
for structural relaxation and electronic structure calcu¬ 


lations. 

In our calculations the substrate Ap © 3 ( 0001 ) was 
modeled by a slab consisting of 6 AI 2 O 3 layers (12 A1 
and six oxygen layers), with the two surfaces terminated 
by Al-I since this termination was found to be more stable 
than the Al-II and the 0-terminated surface^^. Because 
the spacing between the top A 1 layer and the next O layer 
is significantly reduced upon structural relaxatiorP^, the 
slab was firstly fully relaxed. Then a (VTS x V^) sil¬ 
icene layer is placed on each side of a (3x3) supercell 
of Ap © 3 ( 0001 ), which gives rise to a small lattice mis¬ 
match (^2.8%). The silicene/Ap©3(0001) slab is sepa¬ 
rated from its periodic images by ^20 Avacuum regions. 
For the structural relaxation of silicene/Ap©3(0001) the 
Si atoms and the top A1 and © atoms were allowed to re¬ 
lax, with a threshold of 0.001 eV/A for the residual force 
on each atom, while the positions of the other atoms were 
frozen. 

Three types of configurations are considered for 
silicene/Ap©3(0001) as shown in Fig. which are de¬ 
noted by SI, S2, and S3, respectively. S2 (S3) has a Si 
atom at a top site over A1 (©), whereas in SI none of 
Si atoms is exactly sitting on a top site. SI is energeti¬ 
cally more favorable than S2 and S3 by 0.023 eV/Si and 
0.024 eV/Si, respectively, ©ur calculations indicate that 
silicene - substrate distance is increased by about 0.28 
A when vdW correction is taken into account. The re¬ 
laxed structures of silicene are buckled similarly to ideal 
low-buckled silicene, but with significant deviations from 
the alternating buckling pattern of ideal silicene caused 
by interactions with the substrate. For SI the large 
downward bucklings are about 0.9 A, which reduces the 
layer distance down to about 1.4 A. Together with the 
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FIG. 3. Electronic bands of silicene. (a)-(c) /c-projected bands in the (1x1) BZ for isolated silicene in the SI, S2, and S3 
structures, respectively; (e)-(j) the corresponding bands for supported silicene on Al-terminated Al203(0001). (d) and (h) 
/c-projected bands for isolated and supported silicene on 0-terminated Al203(0001), respectively. Bands of ideal(lxl) silicene 
are shown as solid white lines overlaid for comparison. The inset in (a) shows bands about K in the range ±0.04 A“^. The 
Fermi level is set to be zero. 


low-buckled Si atoms, this gives rise to a thickness of 
about 1.5 Afor the silicene. Si-Si bond lengths for the 
relaxed structures are between 2.30 and 2.42 A, com¬ 
parable to those of ideal silicene (2.30 A). The bind¬ 
ing energy of 0.52 eV/Si is comparable to 0.698 eV/Si 
for silicene/Ag(lll)PIlbut much higher than that for sil- 
icene/graphene (^ 0.1 eV/Si), favoring stabilizing the 
monolayer structure. 

The coexistence of low- and high-buckling in the sup¬ 
ported silicene results in a reconstructed structure which 
no longer has the (1x1) silicene periodicity, but still 
maintains the basic honeycomb structure. Therefore, the 
bands can still be unfolded into the (1x1) silicene BZ by 
projecting the supercell wave funct ions onto the corre¬ 
sponding k of the (1x1) silicene celP^^^, i.e., the struc¬ 
tural distortion/reconstruction serves as a perturbation 
to the (1x1) silicene bands. 

Figures. |^a)-(c) show the /^-projected bands for iso¬ 
lated silicene corresponding to the three configurations 
shown in Figs. |^b)-(d), respectively, with the bands of 
the ideal silicene shown for comparison. The isolated (re¬ 
constructed) silicene bands basically preserve the main 
feature of the ideal structure, although there are distinct 
differences. The most prominent one is the opening of a 
gap at the Dirac point at K due to the symmetry break¬ 
ing caused by the structural reconstruction. The size of 
the gap varies for the different structures, and is about 
0.1 eV for the lowest energy structure SI. Moreover, the 
reconstruction gives rise to minigaps in the linear disper¬ 
sion around K, and slightly renormalizes the band widths 
compared to those of the ideal structure, c.f., Fig. |^a). 

Although the bands for silicene on Al-terminated 
Al 203 ( 0001 ), Figs. 3(e)-3(g), are basically preserved, 


there are noticeable changes around the Fermi energy: 
The gap at K for the lowest energy SI structure is en¬ 
hanced to 0.44 eV; additional minigaps are introduced 
(due to the different substrate periodicity); and (com¬ 
paring Figs, ^a) and (e)) the substrate induces a shift of 
the silicene bands of about 0.5 eV relative to the Fermi 
level. In contrast, a comparison of Figs, [^d) and (h) 
shows that silicene bands are significantly disturbed by 
0-terminated AI 2 ©3(0001). In particular, the bands near 
the Fermi level lose features of gapped Dirac cone near 
K as obtained for the isolated silicene. 

The band structure of silicene on AI 2 ©3(0001) 
shows also di stinct differences from metal- 
supported silicen^^®^^ and silicene-based vdW 
heterostructure^^^^^^®^. In the former case, the silicene 
bands are strongly hybridized with those of the substrate 
so that the bands at the K point are significantly de¬ 
stroyed, no longer having Dirac-electron character, while 
the linear dispersions are well preserved in silicene-based 
vdW heterostructures. Moreover, the binding between 
silicene and Al-terminated Al 2 © 3 ( 0001 ), together with 
the structural and electronic properties of the sup¬ 
ported silicene, support that silicene can be properly 
stabilized by a Al 2©3 encapsulation layeil^, although 
the encapsulation layer is likely amorphous in the 
experiment. 

To gain further insights into the interaction between 
the states of silicene and Al2©3(0001), Fig. |^a) shows 
the supercell BZs of silicene/Al2©3(0001), and the BZs 
of the ideal silicene and the substrate. Note that the 
2 / 3 KAI 2 O 3 (KAI 2 O 3 ) coincides with Ksi (KgJ by a trans¬ 
lation of a reciprocal lattice vector of (1x1) Al 2 © 3 ( 0001 ), 
Bi (B 2 ). When silicene is overlaid onto Al2©3(0001), the 
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O-terminated isolated isolated Al-terminated 

Al203(0001) silicene silicene Al203(0001) 


FIG. 4. Interaction between silicene and Al203(0001). (a) 
BZs of silicene/Al203(0001), (1x1) silicene and Al203(0001) 
shown in green, blue and red, respectively. Bi and B2 repre¬ 
sent the reciprocal lattice vectors of (1x1) Al203(0001). (b), 
(c) Band alignment of the isolated silicene and O-terminated 
and Al-terminated (0001) surfaces of AI2O3, respectively. In 
the case of Al-terminated AI203(0001), the isolated systems 
are obtained from configuration SI. The calculated work func¬ 
tions and the gaps of Al203(0001) are shown. VBM and CBM 
represent the valence band maximum and conduction band 
minimum, respectively. 


wavefunctions of silicene 4^si (^Si) interact with those 
of Al203(0001) ^'AbOs (I -f^AbOs+Bi). Figs. I^b) and 
(c) shows the band alignment of the two constituents 
before interacting. In the case of silicene/O-terminated- 
AI 2 03 ( 0001 ) the gapped Dirac states of silicene are close 
to the conduction band of the substrate. This along with 
the strong interaction between them (E5 is larger than 1.0 
eV) lead to strong modifications over the silicene bands. 
Whereas, for silicene/Al-terminated-Al 203 ( 0001 ) (con¬ 
figuration SI) the Dirac states of silicene are in the gap 
of the substrate, far from the substrate states. Therefore 
despite the strong binding between the two constituents, 
there is little direct bonding of between the silicene Dirac 
states and the AI 2 O 3 states. By projecting the wavefunc¬ 
tions over atomic orbitals, we find that Si-pz orbitals con¬ 
tribute over 75% to the VBM and CBM. Whereas, the 
substrate just has little contributions. 

Likewise, a DFT calculation shows that the work func¬ 
tion of an ideal low-buckled germanene is about 4.66 



FIG. 5. DFT calculations of germanene on Au(lll) and 
AI2 03(0001). (a) Relaxed structure of (3x3) germanene on 
(VT^ X a/19) Au(lll). (b) and (c) /c-projected bands for 
isolated and supported germanene, respectively, (d) - (f) 
Gorresponding results for (2x2) germanene on (\/3 x \/3) 
AI203(0001). Buckled Ge atoms are shown in green balls. 


eV, which will place the Dirac point in the gap of Al- 
terminated AI 2 03 ( 0001 ) according to the band align¬ 
ment and will gives rise to an energy separation of about 
1.93 (2.91) eV between the Dirac point and the VBM 
(CBM) of the substrate, respectively. Therefore, one may 
expect that the Dirac states are likely preserved when 
the germanene monolayer is overlaid onto Al-terminated 
AI 2 03 ( 0001 ). DFT calculations were then performed for 
germanene on Al-terminated AI 2 03 ( 0001 ), of which the 
/^-projected bands are shown in Fig. The band struc¬ 
tures for germanene on Au(lll) are also shown for com¬ 
parison. The relaxed structure of germanene on Au(lll) 
is pretty much similar to the planar hexagonal structure. 
Therefore, the /^-projected bands for the isolated ger¬ 
manene are in great similarity to those for the planar 
germanen^^. However, its electronic bands experience 
dramatic changes when germanene is placed on Au(lll). 
From Fig.j^c) one can see that germanene-derived bands 
in the energy window of -6.0 eV to 3.0 eV are strongly de¬ 
stroyed. In particular, near the Fermi level the linear dis¬ 
persions at the K point become indistinguishable. This 
is due to the strong hybridization between the substrate 
bands and the Dirac states. However, germanene on 
AI 2 O 3 exhibits distinct differences. Fig. [^d) shows that 
the relaxed structure of supported germanene basically 
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FIG. 6. Band structure with SOC about the K point for 
germanene/AB03(0001). The bands without SOC (in white) 
are overlaid for comparison. 


remains the global profile of the ideal low-buckled struc¬ 
ture. Therefore it is not surprising that the /^-projected 
band structure for the isolated germanene bears a great 
similarity to that for the ideal germanene, except for 
minigaps near ± 1.0 eV (Fig. [^e)). Upon the presence 
of Al-terminated AI 2 03 ( 0001 ), the linear dispersions near 
K are preserved well, except for the minigap near 1.0 eV 
and the gap at K are slightly enhanced. 

It is demonstrated that external ele ctric field can be 
used to tune the energy gap of silicend^^^. Therefore, 
in addition to band hybridization, electric field induced 
by the substrate may be responsible for the changes in 
the band gaps of silicene and germanene upon supported. 
It should be mentioned that the energy gaps of our sys¬ 
tems obtained using the standard DFT may be under¬ 
estimated. More accurate descriptions can be achieved 
by using hybrid functionals and GW methods. However, 
such calculations are not feasible for our systems with a 
larger number of atoms in the unit cell. Nevertheless, 
our results indicate that the Al-terminated AI 203 ( 0001 ) 
not only stabilizes the monolayer structure of silicene and 
germanene, but also tunes their energy gap for potential 
applications in electronic devices. 

Spin-orbit coupling (SOC) is important for the pre¬ 
dicted quantum spin Hall effect in the ideal low-buckled 
silicene and germanene, which induces a gap opening at 
the K point. The SOC gaps are about 1.55 meV and 
23.9 meV for silicene and germanene, respectively^. Our 
calculations without including SOC show that gap open¬ 
ing occurs in the presence of the substrate (see Fig. 
and Fig. [^, which is due to the symmetry breaking in 
the sublattice. SOC calculations were further performed 
to investigate the effect of SOC on the band structure of 
germanene/Al 203 ( 0001 ). Fig.j^shows that SOC induces 
splittings in both the valence and the conduction bands, 
which are about 25 meV at the K point. Likewise, our 
calculations (plot not shown) for the case of silicene find 
that SOC splittings at K are about 1.5 meV. 

Our results can shed light on the degradation of car¬ 



FIG. 7. /c-projected bands for isolated silicene in the struc¬ 
ture calculated for (a) (Vl2 x \/T2) and (b) (a/TS x \/T3) 
Ag(lll). Bands for ideal silicene are shown as white lines for 
comparison. 


rier mobility observed for silicene-field effect transistors, 
where silicenes directly contact with AI 2 O 3 (although the 
substrate may be amorphous), that is, the measured elec¬ 
tron mobility was an order of magnitude lower than pre¬ 
dicted for the ideal silicen^^. An important issue is that 
the silicenes investigated in the present study differ from 
the experimental ones. In Tao’s experiments, silicene is 
either in (V^ x vT 2 ) or (V^ x vTS) with respect to 
the ( 1 x 1 ) Ag(lll) (which is (V^ x ^/ 7 ) with respect 
to the ( 1 x 1 ) silicene). Fig. shows /^-projected bands 
for isolated silicene in the structures corresponding to 
(VT2 x VT2) and {y/l3xy/l3) Ag(lll). The folded bands 
for the se str uctures agree well those reported by previous 
studie^^^J^. Although Dirac-like bands at K are observed 
in the folded band structure for the isolated silicene on 
(VTSx \/l3) Ag(lll), /c-projected bands indicate that sil¬ 
icene bands undergo considerable changes upon the struc¬ 
tural distortion. The effects of the structural reconstruc¬ 
tion are much more pronounced for the supported silicene 
on (\/T2 X \/l2) Ag(lll). For (\/7 x \/7) silicene placed 
onto AI 2 03 ( 0001 ), a supercell of (\/^ x \/^) is required 
to obtain a small lattice mismatch of ^ 2 %, and a corre¬ 
sponding substrate (V^ x vTS) supercell. While such 
calculations might be possible, we have not attempted 
them; nevertheless, based on the above results one may 
expect similar effects of the substrate on the (^/7 x ^/7) 
silicene as seen here. Moreover, our calculations are con¬ 
sistent with experiments that a gap is formed in the band 
structure of silicene on AI 2 O 3 (Fig. [^. 

The electron mobility is determined by the Fermi ve¬ 
locity and relaxation times of electrons due to various 
scattering processes, i.e., fi = Fig. shows that 

the structural reconstruction renormalizes the electronic 
bands and therefore renormalizes the Fermi velocity as 
well. By a comparison of bands for the ideal silicene 
and the isolated one (Fig. [^a)), a decrease in the Fermi 
velocity for the supported silicene is found. Moreover, 
the substrate interactions and the silicene reconstruc¬ 
tion induce minigaps in the band structure of silicene, 
which also significantly reduce the Fermi velocity near 
the minigaps. The structural reconstructions also mod¬ 
ify the phonon spectrum of silicene and thus the intrinsic 
electron-phonon scattering relaxation time. In addition. 
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surface polar phonon of the substrat e ma y play a role in 
scatterings electrons in the overlayeJ^^^^, further reduc¬ 
ing the electron mobility. 

In summary, we propose a guideline for exploring suit¬ 
able substrates for silicene and germanene. Our DFT 
calculations find that Al-terminated AI 2 03 ( 0001 ) may be 
a good candidate: although silicene (germanene) experi¬ 
ences structural distortions due to the interaction with 
the substrate, the main profile of the ideal low-buckled 
structure is maintained. Unlike the metal-supported 
monolayers, where Dirac electrons are absent, a gapped 
Dirac cone is formed at the K point in silicene as well 
as germanene on Al-terminated AI 203 ( 0001 ). However, 
AI 2 O 3 as a substrate has substantial effects on the elec¬ 
tronic bands of silicene: not only inducing a band gap at 
the K point, but also giving rise to minigaps in the band 


structure, which will lead to Fermi velocity renormaliza¬ 
tion and thus affect the electronic transport properties. 
Our results not only can shed light on recent transport 
experiments on silicene-based transistors, but also may 
aid in the design of silicene/germanene-based electronic 
devices. 

Note added in proof: Since submission of this pa¬ 
per, germanene was grown on defected MoS2, where ger¬ 
manene islands preferentially nucleate at the defect s.l^ 
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